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FOREWORD 


This  is  the  third  and  last  quarterly  report  prepared  under 
contract  AF  1 9(628)  —  594  for  the  Electronic  Systems  Division, 
United  States  Air  Force  System  Command,  L.  G.  Hanscom 
Field,  Bedford,  Massachusetts.  Previous  quarterly  reports 
described  the  technical  performance  agreed  to  and  planned 
under  the  contract.  Full  description  of  the  computer  program 
ESPOD  developed  under  this  contract  appears  in  other  contract 
documentation. 

This  report  was  prepared  by  TRW  Space  Technology  Lab¬ 
oratories,  One  Space  Park,  Redondo  Beach,  California. 
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INTRODUCTION 

This  Third  Quarterly  Report  under  Contract 
AF  19(628)-594  covers  certain  technical  and 
mathematical  considerations  whichare  germane 
to  the  accuracy  of  ESPOD.  Particular  emphasis 
is  given  to  the  numerical  integration  process, 
the  calculation  of  the  vernal  equinox,  and  the  in¬ 
terpretation  of  the  earth  potential  model. 
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1 .  NUMERICAL  EXPERIMENTS  WITH  ESPOD 
INTEGRATION  SUBROUTINE 


This  report  presents  the  results  of  some  numerical  experiments 
conducted  with  the  ESPOD  integration  subroutine. 

The  experiment  consisted  of  integrating  the  equations  of  motion 


(i) 


for  long  periods  of  time  and  then  comparing  the  results  to  reference  orbits. 
The  reference  orbits  were  calculated  by  means  of  analytic  formulas  for  the 
two-body  problem.  In  the  case  of  circular  orbits  this  was  merely  the  ap¬ 
propriate  use  of  a  completely  double  precision  sine-cosine  subroutine. 

For  other  orbits,  Kepler's  equation  was  solved  at  each  compared  time  point. 

A  comparison  between  the  ESPOD  integration  subroutine  and  another 

subroutine  which  is  currently  in  use  in  other  tracking  programs  has  been 

th 

included.  The  ESPOD  integration  subroutine  is  a  10  order  Cowell  scheme 

using  single  precision  arithmetic.  This  routine  was  written  especially  for 

the  Philco  2000  computer  and  the  lack  of  multiple  precision  is  partially 

made  up  for  by  the  larger  word  length  (10  decimal  digits)  of  this  machine. 

A  version  of  this  same  subroutine  was  written  for  the  7094  using  partial 

double  precision.  The  subroutine  called  DE6F,  upon  which  the  Philco 

th 

2000  routine  is  modeled,  is  an  8  order  Cowell  scheme  using  partial 
double  precision  arithmetic.  Besides  the  order  of  the  highest  retained 
differences,  the  other  principal  difference  between  these  subroutines  is 
that  DE6F  restarts  when  it  changes  step-size  while  the  Cowell  method  re¬ 
tabulates  the  differences. 

1 .  1  CASE  1 

The  first  orbit  was  a  circular  orbit  with  r  =  1 .029034544  er,  which  is 

7 

2.  15333073  x  10  feet  or  about  115  miles  above  the  earth.  Such  an  orbit 
will  have  a  period  P  of  88.  19557292  minutes.  Equation  (1)  was  integrated 
for  about  1  6,  000  minutes  or  181  orbits  at  a  fixed  step-size  of  2  minutes. 

The  errors  in  the  individual  rectangular  components  were  sampled' during 
the  first  15  orbits  to  establish  their  pattern  of  growth.  This  showed,  as 
theory  predicts,  that  the  errors  exhibit  an  oscillation  whose  amplitude 
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grows  nearly  linearly  with  time.  The  errors  were  again  sampled  during 
the  last  three  orbits.  The  maximum  error  in  either  of  the  in-plane  com¬ 
ponents  was  less  than  478  ft  for  these  final  orbits.  The  error  in  the  in¬ 
plane  angle 

-1  y2 
0  =  tan  — 

?i 

was  sampled  during  the  entire  run.  The  error  in  0  at  the  end  of  the  run 
-5  -4 

was  about  1x10  radians  or  about  6x10  degrees.  The  same  run  was 
made  with  a  step- size  of  4  minutes  and  the  accumulated  error  in  0  was 
found  not  to  exceed  3.6x10  radians  or  about  2.1x10  degrees . 

1.2  CASE  2 

The  second  orbit  for  which  errors  have  been  computed  is  an  ellipse 
of  eccentricity  e  =  0.73130568  aligned  with  its  major  axis  along  the  y^ 
coordinate  line.  Its  distance  from  the  origin  at  perigee  was  1 . 029865500  er 
and  its  period  was  P  =  634.00369063  minutes.  In  this  case,  Equation  (1) 
was  integrated  until  15  orbits  were  completed  or  for  about  8900  minutes. 
Table  1  gives  the  error  6  observed  in  the  rectangular  components  y^  andy^ 
for  various  subroutines  during  the  14th  orbit.  The  error  is  quoted  in  feet. 
It  was  inconvenient  to  obtain  output  at  exactly  the  same  time  points  on  each 
run  because  the  runs  were  made  in  the  variable  step  mode.  This  accounts 
for  the  lack  of  entries  in  some  places  in  the  table.  In  the  case  of  the  ESPOD 
subroutines,  the  step- size  varied  between  one-half  minute  near  perigee  to 
8  minutes  in  the  vicinity  of  apogee.  In  the  case  of  DE6F  the  variation  was 
one-half  to  4  minutes.  Table  2  gives  a  sample  of  the  errors  6  in  the  rec¬ 
tangular  components  (y^  and  y^)  and  true  anomaly  f  during  the  15th  orbit. 

1.3  CASE  3 

The  third  orbit  for  which  errors  have  been  computed  is  an  ellipse  of 
moderate  eccentricity,  e  =  approximately  0.  015.  The  major  axis  was 
aligned  with  the  y^  coordinate  axis.  Perigee  distance  is  6774  km.  Period 
is  90.  7  minutes.  The  orbit  was  integrated  for  300  revolutions,  that  is, 
for  approximately  18.8  days  or  28,000  minutes.  Reference  points  were 
calculated  for  comparison  from  an  STL  analytic  program.  On  the  last 
orbit  the  maximum  error  in  either  of  the  two  in-plane  components  was 
less  than  5.  7  km. 
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Table  1.  Y-Component  Errors  Observed 
During  14th  Orbit 

ESPOD  on  the  DE6F 


PHILCO  2000 

ESPOD 

on  7094 

(parent  routine) 

T 

(min) 

ISYjft 

|6Y2|ft 

ISYjft 

|6Y2|ft 

IsYjft 

|*Y2|ft 

8244. 5 

176 

1481 

16,  405 

121, 608 

8250 

527 

1298 

349 

727 

8254. 5 

696 

1097 

60,  312 

90,  491 

8256 

815 

705 

552 

400 

8280 

766 

387 

544 

195 

8285 

738 

322 

66,  153 

24, 630 

8295 

678 

224 

504 

73 

82, 594 

15,  571 

8309 

603 

128 

55,  656 

6,  834 

8331 

504 

38 

47,  401 

1,  883 

8355 

416 

23 

349 

123 

8379 

341 

61 

8415 

249 

96 

241 

190 

25,  079 

15,  644 

8475 

126 

136 

146 

213 

8595 

98 

163 

4 

202 

8815 

760 

144 

66,  851 

12,  317 

8850 

938 

655 

381 

287 

82,  402 

58,  864 

Table  2.  Y  Component  and  Phase  Angle  Error 
on  the  15th  Orbit 


T 

(min) 

|sYsjft 

6f  (seconds 
of  arc) 

True  Anomaly 
(deg) 

8865 

420.  8 

747.  1 

6.318 

308 

8876 

-18.9 

971.3 

9.  305 

359.72 

8885 

-424. 1 

821.5 

7.473 

44 

9000 

-353.4 

-153. 1 

0.  645 

149 

9162 

-12.  3 

-247.4 

0.  351 

179.86 

9432 

385.  0 

9.  0 

0.  729 

225 

2.  PRECESSION  AND  NUTATION  OF  THE  EARTH 
FOR  USE  IN  ESPOD 


2.  1  TERMINOLOGY 

2.  1.  1  Luni-solar  Precession 

Long-term  effects  of  Sun  and  Moon  acting  on  the  Earth's  equatorial 
bulge  which  produce  a  conical  rotation  of  the  polar  axis  once  every  26,  000 
years,  or  50"/year. 

2.1.2  Planetary  Precession 

Perturbation  of  the  ecliptic  plane  by  the  planets  causing  a  change  in 
orientation  of  the  ecliptic  about  a  line  at  longitude  174  degrees.  The  effect 
on  the  vernal  equinox  is  about  0'112/year  and  on  the  obliquity  of  the  eclip¬ 
tic  is  about  -0'147/year. 

2.1.3  General  Precession 

The  sum  of  the  effects  of  luni- solar  and  planetary  precession  on 
equator  and  equinox. 

2.1.4  Nutation 

A  relatively  short  period  wobble  of  the  Earth's  pole  superimposed  on 
the  precessional  motion  which  results  from  luni- solar  effects  on  the  equa¬ 
torial  bulge.  The  major  short  period  effect  is  of  amplitude  0'i09  and 
period  13.  6  days. 

2.1.5  Mean  Equator  and  Equinox 

The  plane  perpendicular  to  the  pole  of  the  Earth  including  luni- solar 
precession  is  defined  as  the  mean  equator  at  a  given  date.  The  intersection 
of  this  plane  with  the  ecliptic  plane  including  planetary  precession  of  the 
same  date  defines  the  related  mean  equinox.  The  system  as  a  whole  in¬ 
cludes  general  precession  but  not  nutation. 

2.  1. 6  True  Equator  and  Equinox 

The  plane  perpendicular  to  the  pole  of  the  Earth  including  luni- 
solar  precession  and  nutation  is  the  true  equator  at  a  given  date.  The 


Note  that  the  true  equator  by  definition  corresponds  to  the  actual  equator 
of  the  Earth,  i.  e.,  that  plane  passing  through  the  center  of  mass,  perpen¬ 
dicular  to  the  axis  of  rotation. 
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intersection  of  this  plane  with  the  ecliptic  plane  including  planetary  preces¬ 
sion  defines  the  true  equinox.  The  system  as  a  whole  includes  general 
precession  and  nutation. 

2.  1.7  Of  Date  Versus  of  Date  of  Epoch 

The  above  defined  true  and  mean  systems  may  be  either  rotating  or 
inertial.  If  the  system  is  allowed  to  rotate  with  time  to  include  changes  in 
precession  (and  nutation  for  true  system)  it  is  termed  "of  date."  If  the  sys- 
tern  is  held  fixed  at  conditions  for  some  arbitrary  epoch  (e.  g.,  1950.  0,  0 
of  date  of  launch,  etc.  )  it  is  termed  "of  date  of  epoch,"  or  "of  epoch." 

2.  2  EQUATIONS  AND  ROTATION  MATRICES 

The  fundamental  time  units  and  notations  used  herein  are  as  follows: 


Notation 

Unit 

De  sc  ription 

t 

Ephemeris  day 

Julian  ephemeris  date 

t 

o 

Ephemeris  day 

Julian  ephemeris  date  of  epoch 

*1950.  0 

Ephemeris  day 

Julian  ephemeris  date  of 

1950. 0  =  2433281. 923 

T 

Tropical  century 

(t  -  t  )/36524.220 

T 

o 

Tropical  century 

(to  -  '1950. 0,/36524- 220 

d 

Ephemeris  day 

t  -  2415020. 0  (days  since  1900) 

D 

10,000  ephemeris  days 

(t  -  2415020.  0)/10, 000 

J 

Tropical  year 

(t  -  2415020. 0)/365. 24  (years 
since  1900) 

The  rotation  from  mean  equator  and  equinox  of  epoch  t  to  mean 
equator  and  equinox  of  date  t  is  as  follows: 

r  =  (I  +  P)r 
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where  for  unit  10  : 

Pu  =  -(29,696  +  26To)T2  -  13T3 

P12  =  -P21  =  "(2>  234 -  941  +  1  355Tq)  T  -  676T2  +  22 IT3 
P13  =  "P31  =  -(971.690  -  414To)T  +  207T2  +  96T3 
PZ2  =  -(24,975  +  30To)T2  -  15T3 

P23  =  P32  =  -<10'858  +  2To)T2 

P._  =  -(4721  -  4T  )  T2 
33  '  o' 

The  rotation  from  mean  equator  and  equinox  of  date  to  true  equator 
and  equinox  of  date  is  as  follows 

r  =  (I  +  N)  r 
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where 


( 


N  = 


11 


N11 

N12 

N13\ 

N21 

N22 

N23 

N31 

N32 

N337 

=  0 

"N21 

=  -  AiJj 

COS  C 

:  -N31 

=  -Aiji 

sin  c 

N22  =  0 


N  =  -N  =  -Ae 
23  32  e 


N33  =  ° 


A\|i  =  nutation  in  longitude  of  date 
e  =  obliquity  of  ecliptic  of  date 
Ae  =  nutation  in  obliquity  of  date 

The  nutations  in  longitude  and  obliquity  are  computed  from  sums  of 
sine  and  cosine  terms  of  the  following  form 


Adi  =  Y  S.  sin  (a J  +  bj  '  +  c.F  +  d.D  +  e.fi) 

.1  1  1  1  1  1 
1 

Ae  =  Y  C.  cos  (aj  +  bj  1  +  c.F  +  d.D  +  e.fi  ) 

.1  1  1  1  1  1 
1 

where  the  arguments  (i  ,  l  1 ,  F,  D,  n)  are  functions  of  time  related  to  the 
motions  of  Moon  and  Sun 
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i  =  296?104608  +  1 3?0649924465d  +  0.  0006890D2  +  0?000000295D3 
i  '  =  358?47 5833  +  0.  9856002669d  +  0. 0000112D2  -  0 . 000000068D3 

F  =  11?250889  +  13?2293504490d  -  0.0002407D2  -  0.  000000007D3 

D  =  350?7 37486  +  12?190749194d  -  0?0001076D2  +  0. 000000039D3 

Q  =  259?183275  -  0. 0529539222d  +  0?0001557D2  +  0. 000000046D3 

Note  the  ambiguity  in  notation  where  D  is  both  a  unit  of  time  and  an 
angular  argument.  The  coefficients  S^,  CT  ,  are  grouped  according  to  their 
magnitude  and  period  in  Table  3  with  •  •  ,  e^,  the  integers  multiplying  the 
above  arguments.  Terms  of  magnitude  —01*006  have  been  omitted,  but  may 
be  found  in  Reference  2.5.  1.  Accuracies  have  been  rounded  off  from 
01*0001  (Reference  2.5.  1)  to  01*01. 

For  an  accuracy  reduced  from  Ol'Ol  to  01*1,  the  table  may  be  sum¬ 
marized  by  the  following  expressions 

Aij j  =  -  17 2  sin  ft  +  0'.'2  sin  2ft  -  1 1*  3  sin  (2F  -  2D  +  2ft  ) 

+  0"  1  sin  f  *  -  0 1* 2  sin  (2F  +  2ft  )  +  01*  1  sin  H 

Ae  =  9 V 2  cos  ft  -  0*.’  1  cos  2ft  +  01*6  cos  (2F  -  2D  +  2ft  ) 

+  0 1*  1  cos  ( 2F  +  2ft) 

The  value  for  €  required  to  obtain  sin  €  and  cos  €  multipliers  for 
A^  in  the  N  matrix  is  time  variant.  However,  for  a  period  from  1964  to 
1976  the  change  in  €  is  only  6*.*41  and  for  the  purposes  of  the  computation 
of  the  A \\j  cos  €  and  A^  sin  €  terms  in  the  N  matrix  may  be  taken  as 

€  =  23026'35"  sin  €=  0.39784  cos  €=0.91745 

with  only  a  possible  error  of  ±1  in  the  fifth  decimal  place  for  1964  to  1976. 
For  greater  accuracy,  consult  Table  2.3  of  Reference  2.5.  1,  which  gives 
ecliptic  elements  through  1980. 
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Table  3.  Coefficients  of  Nutation  in  Longitude  and  Obliquity 
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The  rotation  from  mean  equator  and  equinox  of  epoch  to  true  equator 
and  equinox  of  date  is  the  product  of  the  nutation  and  precession  rotations. 
For  application  in  ESPOD  (where  t  -  t^  ^  30  days)  the  second  order  term, 
NP,  is  negligible  and 

£  =  (I  +  N)(I  +  P) 

=  (I  +  N  +  P)  r 


Second  order  terms  in  general  will  not  exceed  an  order  of  2  x  10  T,  where 
T  =  (t  -  t^)/36524.  Hence,  they  may  become  significant  in  the  eight  digit 
after  six  months  and  in  the  seventh  after  five  years. 


The  limiting  magnitudes  of  the  nutation,  30-day  precession,  and  30- 
day  nutation  rotation  matrices  are  presented  below.  The  notation  dN  indi¬ 
cates  the  rotation  from  true  of  epoch  to  true  of  date  (differential  nutation), 
as  opposed  to  N  for  mean  of  epoch  to  true  of  date  (total  nutation). 


/ 

0 

±89 

±39 

N 

T89 

0 

(±)4  9 

\ 

^  =f39 

(=F)49 

0 

(  ° 

-18 

-8 

P  -= 

18 

0 

0 

(30  day) 

8 

0 

0 

(° 

±6 

±3 

dN  -< 

=f6 

0 

(±)3 

(30  day) 

W3 

<=f)3 

0 

\ 


unit  10 


-6 


\ 


The  N32  =  -N23  and  dN32  =  -dN23 
from  the  corresponding  N^,  etc., 

closed  in  parentheses. 


terms  may  be  of  opposite  sign 
terms;  hence,  the  sign  was  en- 
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2.3  COORDINATE  SYSTEMS  USED  IN  ESPOD 


2.  3.  1  Observations 

Radar  observations  are  by  virtue  of  their  connection  to  the  actual 
position  of  the  Earth  in  true  equator  equinox  of  date  system.  Baker-Nunn 
angular  observations  are  available  to  ESPOD  in  varying  coordinates  de¬ 
pending  on  whether  they  are  precision  or  field  reduced  and  whether  or  not 
they  are  preprocessed  by  other  programs  in  the  SPADATS  system  before 
entering  ESPOD.  The  applicable  logic  is  given  in  the  introduction  to 
Section  7.  1  of  ESD- TDR-64- 395  and  the  updating  within  ESPOD  is  de¬ 
scribed  under  the  ADJUST  subroutine  in  the  same  document. 

2.3.2  Integration 

The  integration  is  carried  on  using  the  time  equator  and  equinox  of 
0  •  0  on  the  day  of  epoch.  The  right  ascension  of  Greenwich  at  epoch  is 

Qgo  =  43206  I23925-  836  +  1  •  84542 J  +  9.29  x  10"6JZ  +  n  +  2tt(j  -  [j 

r  _lo  '  t1900.  0 
J  "  — 3^.TB — 


greatest  integer  in  J 


The  expression  n  changes  a  to  true  equator  and  equinox  of  0  •  0  of  the 

day  of  epoch.  Both  the  rotation  of  the  orbit  to  Earth-fixed  coordinates  for 
representation,  and  the  computation  of  station  longitude  are  dependent  on 
a  for  determining  the  position  of  the  Earth  at  a  given  time. 


The  computation  for  a  includes  the  nutation  correction  which 

g° 

changes  the  station  right  ascensions  and  places  the  integration  in  true  of 
0  •  0  of  day  of  epoch  (see  subroutine  TINIT).  The  precision  and  field- 

reduced  angular  observations  are  also  corrected  for  nutation.  The  correc- 
tion  is  made  using  A^,  Ae  ,  at  0  •  0  of  day  of  epoch,  to  agree  with  the 

integration  (see  subroutine  ADJUST).  The  adjusted  values  of  the  precision 
reduced  angular  observations  will  vary  from  true  of  epoch  by  a  negligible 
P(t  -  t^),  the  precession  from  0  •  0  of  day  of  epoch  to  the  time  of  observa¬ 

tion  because  of  their  original  reference  to  mean  of  date. 
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2.  3.  3  Summary  of  Coordinate  Systems 

The  coordinate  systems  after  revision  will  be  as  follows: 


Radar  observations 

Field-reduced  observations 
Baker-Nunn 

Precision  reduced  Baker-Nunn 
observations 

Integration 

Station  longitude 

Station  latitude 

2.3.4  Bias  Limits 


true  of  date 
h 

true  of  0  *0  day  of  epoch 

Vi 

true  of  0  *0  day  of  epoch  +  P(t  -  t^) 

true  of  epoch 
true  of  epoch 
true  of  date  =  true 


In  observed  minus  computed  residuals 


Type  of 
Observation 


Bias 


P(t  -  t  )  +  dN 
o 


Radar 

Field  reduced  angular 
Precision  angular 

2.4  APPROXIMATE  NUTATION  FORMULAS 


None 


p(t  .  toi 


Limiting 

30 -Day  Magnitude 

0.15  km 
0 

0.11  km 


Nutation  in  longitude,  and  obliquity  at  epoch,  are  obtained  to  within 
01' 5  from 

A^  =  -  1 7 2  sin  Q  -  1 »  3  sin  (2F  -  2D  +  2ft) 

Ae  =  9 "2  cos  ft  +  01*6  cos  (2F  -  2D  +  2ft ) 

£2  =  259?18  -  19?3414  J 
2F  -  2D  +  2ft  =  -l60?6l  +  7 1 9?9957  J 


The  correction  n  to  a  is  then 

g° 


n  =  cos  c  Av|j 
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and  the  corrections  for  nutation  to  the  precision  and  field- reduced  Baker- 
Nunn  observations  are 

Aa  =  (cos  €  +  sin  e  sin  a  tan  6)A^i  -  cos  a  tan  6Ae 
A6  =  sin  €  cos  aAijj  +  sin  aAt 
2.5  REFERENCES 

2.5.  1  Explanatory  Supplement,  American  Ephemeris  and  Nautical  Almanac, 

1961. 

2.5.2  ESD-TDR-64- 395  ESPOD  Mathematical  and  Subroutine  Description. 


\ 
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3.  GEOPOTENTIAL  MODEL  ACCURACIES 


3.  1  ACCURACIES  OF  SPADATS  MODELS 

Four  coherent  models  of  the  earth  potential  field  will  be  provided 
with  ESPOD.  These  models  are  summarized  in  Table  4. 

3.  1.  1  Model  1 

Model  1  will  give  an  accuracy  within  approximately  0.  5  km  of  a 
period  of  24  hours.  The  major  error  is  caused  by  neglecting  the  second 
sectorial  harmonic  (with  coefficient  ^  at  t^ie  phase  angle  X ^  .  The 

effect  is  sinusoidal  with  a  period  of  12  hours. 

3.1.2  Model  2 

Model  2  will  give  approximately  the  same  accuracy  as  Model  1  for 
the  same  reason. 

3.  1.  3  Model  4 

Model  4  is  estimated  to  provide  an  accuracy  of  approximately  0.  2  km, 
with  uncertainties  in  the  sectorial  and  tesseral  harmonics  making  further 
accuracy  impossible. 

3.2  REFERENCES 

3.2.1  R.R.  Newton,  Results  from  Transit  Satellite  Tracking,  reported 
in  A.H.  Cook,  Space  Science  Reviews,  2,  1963. 

This  report  indicated  an  amplitude  0.  5  km  for  the  ^  perturbation. 

3.2.2  R.R.  Newton,  Science ,  15,  803-808. 

A  lower  error  limit  of  0.  10  km  is  estimated  possible  by  eventual 
inclusion  of  50  harmonic  coefficients  in  a  solution  with  ample  and  well 
distributed  data. 

3.  2.  3  Conclusions  referenced  in  Paragraph  3.  1.  3  are  an  estimate  by 
M.P.  Francis  from  probable  errors  associated  with  Model  4  by  Kozai, 

SAO  Special  Report  No.  72,  9  August  1961. 
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Table  4.  Earth  Potential  Model  Options 
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3.  3  TESTS  FOR  POTENTIAL  HARMONICS 


Tests  have  been  run  to  determine  the  effect  of  omitting  higher  order 
harmonics  from  the  earth's  gravitational  field.  Although  the  tests  were 
not  run  for  the  geopotential  coefficients  provided  with  ESPOD,  they  are 
easily  converted  to  that  reference  and  provide  in  any  case  an  estimate  of 
the  relative  influence  of  various  geopotential  terms. 

The  test  was  performed  by  taking  initial  conditions  for  a  100  nautical 
mile  altitude  circular  orbit  over  a  (hypothetical)  spherical  earth  and  cal¬ 
culating  for  that  orbit  the  partial  derivatives  of  the  cartesian  coordinates 
with  respect  to  the  coefficients  of  the  geopotential  harmonics.  The  partial 
derivatives  at  the  points  of  1/4  orbit,  1/2  orbit,  3/4  orbit  and  one  full 
orbit  were  root-sum  squared  to  obtain  the  total  effect  of  the  coefficient  on 
each  position  and  velocity.  These  parial  derivatives  were  then  multiplied 
by  the  values  of  the  coefficients  and  uncertainties  tabulated  in  Table  5. 
These  values  were  derived  from  W.M.  Kaula's  report,  "NASA  Technical 
Note,"  D- 1848,  June  1963. 

Two  cases  were  run  with  the  defined  orbit,  one  starting  at  the  north 
pole  on  a  polar  orbit  and  one  starting  at  the  equator  on  a  polar  orbit. 

Tables  6  and  7  show  the  results  of  this  study. 


> 
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Table  5 


I 


.  Geopotential  Coefficients  from  NASA  Technical  Note 


D- 

1848, 

June  1963 

n 

m 

J  X  106 

n,  m 

X 

n,  m 

2 

1082. 30  ±0.01 

3 

-2. 59  ±0. 03 

4 

-1. 65  ±0.  30 

5 

-0. 10  ±0. 03 

6 

0.  36  ±0. 07 

7 

-0. 39  ±0. 04 

2 

2 

1. 62  ±0.  15 

-2 1?4  ±2?8 

3 

3 

0.11  ±0.06 

37?6  ±11?0 

4 

4 

0.013  ±0. 006 

28?4  ±6.  5° 

3 

1 

1. 92  ±0.  23 

-3?6  ±6?5 

3 

2 

0. 12  ±0. 09 

±28?4 

4 

1 

0. 48  ±0. 14 

±17?9 

4 

2 

0.072  +0. 042 

47?7  ±16?9 

4 

3 

0. 031  ±0. 012 

5?9  ±7?0 
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Table  6.  Effect  of  Gravitational  Harmonics  and  Their  Uncertainties 
on  a  Polar  Orbit  Starting  at  North  Pole 


At  One  Fourth  Revolution 

At  One  Half  Revolution 

Position,  feet 

Velocity,  ft/sec 

Position,  feet 

Velocity  ft/ sec 

J2 

147,  330  ±1.4 

164. 71  ±0.002 

272,  900  ±2.  5 

322. 42  ±0. 003 

J3 

261.  9  ±3.  0 

0.  30  ±0. 003 

636.  9  ±7.  4 

0. 65  ±0. 008 

25.  8  ±4.  7 

0.01  ±0.  002 

416. 2  ±75. 7 

0. 49  ±0.  09 

4.  6  ±1.4 

0. 002  ±0. 0005 

24. 3  ±7.  3 

0. 02  ±0. 007 

*6 

17. 5  ±3.4 

0.01  ±0.  002 

87. 8  ±17. 1 

0.  10  ±0. 02 

J7 

22. 1  ±2. 3 

0. 03  ±0. 003 

91.8  ±9.  4 

0. 09  ±0. 009 

/(2) 

2 

882.  8  ±81.7 

1.  12  ±0.  10 

455. 5  ±42. 2 

0. 52  ±0. 05 

j3n) 

694.  8  ±83.  2 

0. 74  ±0.  09 

389.0  ±46. 6 

0. 46  ±0.06 

J3(2) 

37. 5  ±28.  1 

0.01  ±0.01 

52. 6  ±39. 5 

0. 04  ±0. 03 

j3(3) 

66.8  ±36 . 4 

0. 08  ±0.  04 

184. 2  ±100. 5 

0. 20  ±0. 11 

J4(l) 

83. 4  ±24. 3 

0.  03  ±0. 008 

52. 3  ±15. 3 

0.  09  ±0.  03 

j4(2) 

66. 2  ±38. 6 

0. 09  ±0.  05 

32. 3  ±18. 9 

0. 03  ±0. 02 

J4(3) 

74. 9  ±29. 0 

0.  12  ±0.  05 

89.  9  ±34.  8 

0.  08  ±0.03 

J4(4) 

60. 9  ±28.  1 

0. 12  ±0. 06 

192. 3  ±88. 8 

0. 22  ±0. 10 

At  Three -Fourths  Revolution 

At  One  Revolution 

Position,  feet 

Velocity,  ft/sec 

Position,  feet 

Velocity,  ft/sec 

J2 

481, 030  ±4.  4 

515. 94  ±0.005 

520, 790  ±4. 8 

616. 56  ±0. 006 

^3 

867. 6  ±10. 0 

0.  92  ±0.01 

731. 5  ±8. 5 

1.01  ±0.01 

J4 

728. 3  ±132. 4 

0. 78  ±0.  14 

797. 1  ±144. 9 

0.  94  ±0.  17 

J5 

31.2  ±9.4 

0.  03  ±0.01 

24. 7  ±7.  4 

0.  03  ±0.01 

J6 

153.  3  ±29.  8 

0.  16  ±0.  03 

168.  1  ±32. 7 

0. 20  ±0.  04 

^7 

115.3  ±11.8 

0.  12  ±0.01 

87. 7  ±9. 0 

0.  13  ±0.01 

j2<2> 

817. 2  ±75. 7 

0. 95  ±0.  09 

1055. 8  ±97. 8 

1. 22  ±0.  11 

J3IU 

388. 6  ±46. 6 

0. 32  ±0.04 

90. 9  ±10. 9 

0.  11  ±0.  01 

J3(2) 

36. 9  ±27. 7 

0.06  ±0. 04 

85. 3  ±63. 9 

0. 05  ±0.  04 

j3(3) 

271. 4  ±148. 0 

0. 24  ±0. 13 

141. 2  ±77. 0 

0.  15  ±0.  08 

J4(!) 

106. 6  ±31. 1 

0.0  3  ±0.01 

11.5  ±3.  4 

0. 15  ±0.04 

j4(2) 

68. 9  ±40. 2 

0. 09  ±0.  05 

115.2  ±67.  2 

0.  13  ±0. 07 

J4(3) 

85. 4  ±33. 1 

0. 07  ±0. 03 

17. 9  ±6. 9 

0.  13  ±0.  05 

j4<4> 

354. 9  ±163. 8 

0.41  ±0.19 

474. 3  ±218. 9 

0. 54  ±0.25 
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)  Table  7.  Effect  of  Gravitational  Harmonics  and  Their  Uncertainties 

on  a  Polar  Orbit  Starting  at  Equator 


At  One  Fourth  Revolution 

At  One  Half  Revolution 

Position,  feet 

Velocity,  ft/sec 

Position,  feet 

Velocity,  ft/sec 

J2 

22, 380  ±0.2 

22. 07  ±0. 0002 

56, 141  ±0. 5 

66.  98  ±0. 0006 

'3 

21.7  ±0.3 

0. 05  ±0. 0006 

106.0  ±1.2 

0. 10  ±0. 001 

J4 

29.  9  ±5.  4 

0.09  ±0.02 

229. 3  ±41. 7 

0. 27  ±0.  05 

■^5 

1.4  ±0.4 

0.004  ±0.001 

8. 7  ±2. 6 

0.01  ±0. 003 

J6 

4. 2  ±0.8 

0.005  ±0. 001 

4. 0  ±0.8 

0. 005  ±0. 0009 

^7 

3.1  ±0.3 

0. 007  ±0. 0008 

20. 7  ±2. 1 

0. 02  ±0. 003 

J2(2) 

291.4  ±27.  0 

0.47  ±0.  04 

1,  173. 0  ±108. 6 

1. 44  ±0. 13 

J3U) 

94.0  ±11.3 

.0.07  ±0.  009 

502. 4  ±60. 2 

0. 60  ±0. 07 

J^2) 

35. 1  ±26. 4 

0. 06  ±0.04 

111.2  ±83.  4 

0. 13  ±0. 10 

J3(3) 

88.2  ±48. 1 

0.  11  ±0. 06 

228.7  ±124. 8 

0. 20  ±0. 11 

J4(1) 

7. 9  ±2. 3 

0.04  ±0.01 

71.3  ±20.8 

0. 04  ±0. 01 

J4(2) 

12. 0  ±7.0 

0.04  ±0.  02 

18.  1  ±10. 6 

0.04  ±0. 02 

j4<3> 

67.0  ±25. 9 

0. 10  ±0. 04 

142. 2  ±55.0 

0. 09  ±0. 04 

J4I4> 

85.6  ±39.  5 

0.  11  ±0.05 

262.6  ±121. 2 

0.35  ±0.  16 

At  Three -Fourths  Revolution 

At  One  Revolution 

Position,  feet 

Velocity,  ft/ sec 

Position,  feet 

Velocity,  ft/  sec 

J2 

101, 350  ±0.  9 

105. 41  ±0.001 

102,  630  ±0. 9 

122. 67  ±0.001 

^3 

225.5  ±2. 6 

0. 17  ±0. 002 

121. 6  ±1.4 

0. 15  ±0. 002 

405. 0  ±73.  6 

0. 17  ±0. 03 

421. 9  ±76. 7 

0.  50  ±0.  09 

12.4  ±3.  7 

0. 002  ±0.0007 

5.  6  ±1.7 

0.007  ±0. 002 

*6 

8.2  ±1.6 

0. 009  ±0. 002 

6.  1  ±1.2 

0. 007  ±0. 001 

J7 

42.4  ±4. 3 

0.03  ±0. 003 

22. 5  ±2. 3 

0. 03  ±0. 003 

J2(2) 

2, 135. 2  ±197. 7 

2. 32  ±0.21 

2,  431. 6  ±225. 2 

2.90  ±0.  27 

J^1) 

1065. 5  ±127. 6 

1. 22  ±0.  15 

1248.8  ±149. 6 

1. 25  ±0. 15 

J3(2) 

124. 5  ±93. 3 

0. 11  ±0.08 

43. 2  ±32. 4 

0.06  ±0 .  04 

J3(3> 

192. 7  ±105. 1 

0. 18  ±0. 10 

34. 6  ±18.8 

0.  16  ±0.  09 

V‘> 

10.3  ±3.0 

0. 12  ±0. 04 

132.5  ±38.6 

0. 001  ±0. 0004 

j4<2> 

56. 6  ±33. 0 

0.06  ±0. 04 

87. 9  ±51. 3 

0. 11  ±0.  07 

J4(3) 

70.6  ±27.3 

0. 12  ±0.05 

113.  1  ±43.8 

0.03  ±0.01 

J4(4) 

554.2  ±255.8 

0.61  ±0.28 

700. 9  ±323. 5 

0.82  ±0.  38 
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4.  APPLICATION  AND  THEORY  OF  THE  A  PRIORI  ATA  MATRIX 

(S- MATRIX) 


4.  1  ESPOD  OUTPUT 


ESPOD  will  punch  out  the  A  A  matrix  evaluated  at  epoch  following 

a  differential  correction  by  selecting  this  option  with  JDC  card  column  44 

T 

set  to  1.  ESPOD  will  punch  out  the  A  A  matrix  evaluated  at  an  arbitrary 
time  following  a  differential  correction  by  selecting  this  option  with  JDC 
card  columns  51,  55  and  56  set  to  1,  assuming  a  correct  specification  of 
the  arbitrary  time. 


4.  2  GENERAL  APPLICATION  OF  A  PRIORI  ATA  MATRIX 
T 

The  A  A  matrix  summarizes,  scales,  and  weights  the  contribution 
of  a  set  of  residuals  to  the  differential  correction  for  a  set  of  orbital  ele¬ 
ments  after  an  integration  through  a  set  of  data  over  a  given  time  span. 

T 

The  A  A  matrix  is  specific  to  the  time  at  which  the  orbital  elements  are 

defined.  If  it  is  desired  to  use  the  same  residuals  to  affect  a  subsequent 

correction  to  the  elements,  in  combination  with  a  newly  observed  set  of 

T 

residuals,  the  punched  out  A  A  matrix  permits  this  application.  The 

time  at  which  the  elements  to  be  corrected  will  be  defined  must  be  known 
T 

and  the  A  A  matrix  for  the  earlier  residuals  developed  for  that  time.  The 
T 

A  A  matrix  is  then  provided  to  ESPOD  with  preliminary  data  and  by  set- 

❖ 

ting  JDC  card  column  42  equal  to  1 .  See  Reference  2.5.2,  Section  5-12 

T 

for  further  details  and  restrictions  on  the  use  of  the  a  priori  A  A  matrix. 


4.  3  THEORY  (WITH  EXAMPLE) 
4.  3.  1  Run  1 


Let  elements  E^  define  residuals  Rj  from  observations  Oj  and  also 

T  _ 

define  AjAj  =  Sj.  Together  they  define  a  correction  AEj.  Let  elements 

E2  =  Ei  +  AEj,  that  is,  those  elements  corrected  from  E  \  by  residuals 

T 

R  i  on  Run  1.  Let  AjAj  be  punched  out. 


See  Reference  2.  5.  2,  Section  6.  1.  4,  for  the  definition  of  A  A. 
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4.3.2  Run  2 


Let  elements  E£  define  residuals  R 2  from  observations  O;?  and  also 
T  T 

define  A2A2.  Use  A^Aj  as  an  a  priori  matrix,  which  during  Run  2 

T  T  — 

accumulates  in  the  sum  AjA^  +  A2A2  =  S2*  Elements  E2,  together  with 
R2  and  S2,  define  a  correction  AE2.  Let  elements  E3  =  E2  +  AE2. 

4.3.3  Check  Run  3 

Let  elements  E^  define  residuals  R\  and  R2  from  observations 

T 

Oj  and  C>2  and  also  define  A3A3.  Together  they  define  AEj.  Let  elements 
E3  =  E^  +  AE^ 


4.3.4  Proof 

See  Reference  2.  5.  2,  Section  6,  for  general  development  of  this 
notation. 


Run  1 


£j(E) 

A1 

f2(E) 


Weighted  residual  vector,  data  1 
Partials,  data  1 

Weighted  residual  vector,  data  2 
Partials,  data  2 


A[Al(E2  -  E,)  = 


Run  2 


(A7Al  +  AJa2)(e3  -  E2)  .  Ajf2(E2)  =  aJe2(ei)  -  aJa2(e2  -  EJ 


(aTa, 


+  a^a2) 


E3  -  E1  -  (E2  -  El)i  =  Ak  (El)  -  A2A2  (e,  -  Bj 


KA1  +  a2ta2)  (e,  -  eJ  -  aJa,  (e2  -  eJ  -  a^a2  (e2  -  eJ 

(El)  '  A2A2  (E2  -  El) 


-  ATf 
~  A2i2 
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KA1  +  AIA2)(E3  -  El)  =  AIf2(E2)  +  A[Al(E2  -  El) 

KA1  +  A2A2)(e3  -  Ej  =  AJ«2(e2)  +  aJtJe,) 

Run  3 

(AfAl  +  A2A2)  (E3  -  El)  =  A?£l  (El)  +  AJf2  (E2) 
Thus,  '= 


i 
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